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ABSTRACT 

We describe the discovery of a bright, young Kuiper belt-like debris disk around HD 115600, a ^ 

1.4-1.5 Mq, ^ 15 Myr old member of the Sco-Cen OB Association. Our H-band coronagraphy/integral 
field spectroscopy from the Gemini Planet Imager shows the ring has a (luminosity scaled) semi major 
axis of 22 AU) ^ 48 AU, similar to the current Kuiper belt. The disk appears to have neutral 
scattering dust, is eccentric (e ^ 0 . 1 - 0 . 2 ), and could be sculpted by analogues to the outer solar system 
planets. Spectroscopy of the disk ansae reveal a slightly blue to gray disk color, consistent with major 
Kuiper belt chemical constituents, where water-ice is a very plausible dominant constituent. Besides 
being the first object discovered with the next generation of extreme adaptive optics systems (i.e. 
SCExAO, GPI, SPHERE), HD 115600’s debris ring and planetary system provides a key reference 
point for the early evolution of the solar system, the structure and composition of the Kuiper belt, 
and the interaction between debris disks and planets. 

Subject headings: planetary systems, stars: solar-type, stars: individual: HD 115600 


1. INTRODUCTION 

Situated just beyond Neptune’s orbit (30-50 AU), the 
Kuiper belt is home to thousands of remnants of the 
earliest stages of icy planet formation, is the location 
of numerous dwarf planets such as Pluto, and provides 
keys to understanding the early, unevolved solar system 
(|Jewitt and Luul[l992l) . Cold debris rings around nearby, 
young stars offer a critical refe rence for the Kuiper 
belt’s evolution an d composition (iKenvon and BromlevI 
120081 : iWvattI 1200811 . However, the few such rings im¬ 
aged in scattered light surround stars at stellocentric 
distances farther than the Kuiper belt (e.g. iKalas et al.l 

120051 : ISchneider et al.l[2009t IKrist et 3111201^ , are not lo¬ 

cated in massive OB association resembling th e solar sys¬ 
tem’s likely birth environment (IAdamsll 2010 f) . and lack 
spatially-resolved scattered-light spectroscopy to explore 

composition. __ 

HD 115600 {d = 110.5 pc: Ivan Leeuwenll200^ is a 1.4 
M 0 , F2V/F3V 15 Myr-o\d member of the Lower Cen- 
taurus Crux (LCC) region of the Sco-Cen OB Associ¬ 
ation (|Pecaut et al.l 1201211 0. The star has a large in¬ 
frared excess with a dust mass and fractional luminosity 
(0.05 Mmoon; Lir/L* ^ 1.7 X10“^) comparable to the 
well-known lumin o us, resolved ring aroun d HR 4796A 
(|Chen et al.l 1201 iL 120151 : lJura et all 1199811 . Thus, HD 
115600 is a promising target for high-contrast imaging. 
We report the discovery of a bright debris ring around 
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HD 115600 using integr al field spectroscopy fr om the 
Gemini Planet Imager (|Macintosh et al.l 1201411 . The 
ring is plausibly sculpted by an unseen solar system-like 
planet, surrounds the star at a Kuiper belt-like distance, 
and has a spectrum consistent with major Kuiper belt 
constituents. The HD 115600 planetary system provides 
a valuable reference point for the early evolution and 
composition of the solar system, the structure of the 
early Kuiper belt, and the interaction between debris 
disks and planets. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained integral field spectroscopy of HD 115600 
on 23 March 2014 w ith the Gemini Planet Imager 
(jMacintosh et al.ll2014l l in H band (Ao = 1.6 ^m, R = 44- 
49; 07014 pixeU^) behind the apodized pupil lyot corona- 
graph (rm.n.s k ^ 0712) and in an gular differential imaging 
mode (API: iM^arois et akl 120061 1. Our observations con¬ 
sisted of 58 co-added 49.5 second frames, yielding a total 
integration time of ~ 48 minutes and a field rotation of 
29° (3 A/D at 0725). Ae and Ar lamp observations pro¬ 
vided wavelength calibration. 

We carried out basic data reduction usi ng the Gem¬ 
ini D ata Reduction Pipeline, version 1.2.1 (|Perrin et al.l 
l2014fl . We modifie d the softwa r e, add ing image regis¬ 
tration steps as in iGurrie et al.l (|2011h . Removing the 
middle 20 data cubes, which exhibited a pupil stop mis¬ 
alignment, and those with poorer AO correction reduces 
our total integration time to 28 minutes. 

For poi nt-spread function (PSF) subtraction, we used 
A-LOCI (IGurrie e t al.ll201 ^ l2014cll tu ned to “conserva- 
tive” settings (jThalmann et al.l 120111) better suited for 
recovering extended emission, imposing a rotation gap 
criterion of 6 = 1.5, an optimization area of Na = 
1000, and a singular value decomposition (S VD) cutoff of 
STDiirn = 10“® (see iLafreniere et al.l[2007l : IGurrie et al.l 
l2014bf) . We construct a final data cube from a median- 
combination of PSF-subtracted cubes and wavelength- 
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collapsed using an outlier-resistant mean to yield a band- 
averaged image. Due to the small nu mber of data cubes , 
we did not employ speckle filtering (|Currie et al.ll201§ . 
l2014bH to truncate the set of reference PSFs. For initial 
flux calibration, we measured the average, background- 
subtracted signal of the four satellite spots (transmission 
= 2.035 X 10“^) in an 8-pixel diameter aperture and 
divided the data c ube by a normalized F2V spectrum 
(jPickles et al.llTo^ . 

3. DETECTION OF THE HD 115600 DEBRIS DISK 

Figure [T] (top panels) displays our wavelength- 
collapsed disk image and an individual channel, clearly 
revealing a bright, thick debris ring visible from r ~ 0"25 
to r ^ 0"55 and viewed close to edge-on. The disk ansae 
along the major axis extend from r ^ 0'.'34 to 0"5, or 
37.5 AU to 55 AU, comparable to the current Kuiper 
belt. Both sides of the disk are visible, ruling out strong 
forward-scattering anisotropy. The disk is offset from the 
star position and vi sually resembles a puf fier version of 
the HR 4796A disk (|Schiieider et al.ll2009ll . 

To estimate the disk’s SNR in the collapsed image and 
in each spectral channel, we consider two methods. First, 
for the colla p sed im age we construct a SNR map as in 
iCurrie et al.l (1201111 . replacing each pixel with the sum 
of values enclosed by a 3-pixel aperture (^ 1 FWHM) 
and computing the robust standard deviation at each 
angular separation. For a given region of the disk, this 
procedure considers other disk regions (e.g. those on the 
opposite side) and negative self-subtraction residuals as 
“noise”, which leads to underestimating the significance 
of the disk signal compared to real residual speckle noise. 
Nevertheless, the disk ansae in the collapsed image are 
still detected at SNR ^ 8; the disk’s signal exceeds 2-cr 
exterior to r ^ 0'.'25. 

Second, we compute the ansae SNR in each spectral 
channel, defining the noise likewise as the robust stan¬ 
dard deviation of the summed image in a FWHM-wide 
arc at the same angular separation as the disk but at 
different position angles, at a separation (10-50 pixels 
away) chosen to avoid the disk signal and negative self¬ 
subtraction footprints. The SNR in each spectral channel 
in the disk ansae ranges between 5 and 8; the SNR of the 
collapsed image exceeds 10 in the ansae. 

4. ANALYSIS 

To derive the disk geometry and extract spectra, we 
follow a three-step approach. First, we derive the de¬ 
bris ring’s basic geometry from the wavelength-collapsed 
data cube using ellipse fitting. Second, we use forward¬ 
modeling to fine-tune these properties and calculate 
second-order properties of the disk (e.g. scattering func¬ 
tion). Third, we extract the disk’s spectrum in the bright 
ansa regions, using the best-fit model to the final data 
cube to calibrate the disk throughput in each spectral 
channel. 

4.1. Geometry 

To determine the basic geometry of the disk - in¬ 
clination, position angle, semi-major/minor axes, and 
disk center - we first use the using the IDL mpfitel- 
lipse package to determine an approximate trace of the 
disk, where the pixels are weighted by their SNR. Sec¬ 
ond, we constructed a grid of ellipse parameters around 


the best-fit set determined by mpfitellipse, and calcu¬ 
lated a more robust best-fit value using the “maximum 
merit” procedure, identifying the ellipse parameters that 
maximized the disk sig nal along the trace of the ellipse 
(jThalmann et al.ll2011fl . We repeat this step using dif¬ 
ferent ranges in radii and different cutoffs in SNR for the 
disk trace (e.g. SNR > 1.5, 2; r = 17.5-45 pixels; 20- 
40 pixels) to define best-estimated values and associated 
uncertainties. 

We derive a best-fit position angle and inclination of 
PA = 24° ± 0.5° and i = 79.5° ± 0.5° (1-cr errors). The 
disk semimajor/minor axes are rmajor minor = 0'.'425 ± 
O'.'OlO (- 48 AU ± 1.1 AU), 0"077 ±’0'.'007 (8.5 ± 0.8 
AU). The semimajor axis and visible extent of the disk 
ansae are within the range of radii encompassing the 
current Kuiper belt; the luminosity-scaled (rmajor/VT^) 
semi major axis (extent of the ansae) are ^ 22 AU (17- 
25 AU), comparable to the predicated distances of ma¬ 
jor models of t he early, pre-stirred Ku iper belt,: e.g. 
the N ice model (iLevison et al.l 1200811 or iNesvornv et al.l 
(|2015ll . The projected disk center is offset from the star 
Ax,Ay = 0"018 ± 0"008, 0"029 ± 0"014 (Fig. 1; dia¬ 
monds) and the peak signal of the two ansae differ by ^ 
0.5 pixel in angular separation (~ 1 AU). 

4.2. Disk Forward Modeling 

To infer additional disk properties, we generate a grid 
of sy nthetic scattered light images using the GRaTeR 
code (lAiigereaii et al.l 1199911 , and forward-model these 
synthetic disks using our pipeline to compare the pro¬ 
cessed synthetic disk image with the real disk image in 
each spectral channel, extending to integral field spec¬ 
troscopy met hods that have been applied to broadban d 
imaging data (jEsposito et al.l[2?)T^ iMazover et al.l[2?iT3l . 
We insert a disk model into a sequence of empty data 
cubes with position angles identical to those of our sci¬ 
ence sequence, and convolve the model in each spectral 
channel with the appropriately-sized PSF. We performed 
PSF subtraction on the cubes containing the model disk 
using the same A-LOCI coefficients that were applied to 
the real data. 

Table m describes the model parameter space we ex¬ 
plore and our results. To reduce the dimensionality of 
our forward-modeling, we adopt the inclination and po¬ 
sition angle determined from our ellipse modeling (79.5° 
and 24°, respectively). We defined the argument of peri- 
center to be [90, 270]°, or 90° from the visible disk major 
axis: departures from these values resulted in brightness 
asymmetries between the disk ansae inconsistent with 
the data. We varied other parameters, including the the 
ring center (oo = 47 and 48 AU), the Henyey-Greenstein 
scattering parameter g (0-0.15), the disk scale height at 
the disk center (ksio = 1-5 AU), the density power laws 
describing the decay of ring emission away from the ring 
center (am = 7.5, 10; aout = —5, —7.5), the disk eccen¬ 
tricity (0-0.3), and the offset along the major axis (0 and 
1 AU). While our parameter space search is not exhaus¬ 
tive, values outside these ranges (e.g. g > 0.15) yielded 
processed synthetic disk images strongly discrepant with 
the real data. 

To match the observed brightness in each spectral 
channel, we fix the parameters in the model and scale 
the flux and define the fit of the model to the data us¬ 
ing the real and synthetic collapsed images convolved by 
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the PSF as in our SNR calculations. Our “acce ptably fit¬ 
ting” models are those f ulfilliuK < Xmm + x Ndata 
(|Thalmann et al.l 1201311 . We perform our optimization 
(and define x^) by the visible trace of the disk at sepa¬ 
rations where it is consistently visible at SNR > 2 along 
both sides 0"28 and 0'.'52). 

The best-fit model is an offset, eccentric ring with 
neutral-scattering dust (Figure [T] (bottom panels)). This 
model exhibits some discrepancies with the data, slightly 
over (under) predicting the signal on the southeastern 
(southwestern) side. However, our best-fit model gen¬ 
erally reproduces disk shape, providing a good first in¬ 
vestigation. More sophisticated models fit to higher SNR 
data and exploring more parameter space will better con¬ 
strain disk properties. Additionally, much of our param¬ 
eter space remains degenerate, admitting widely vary¬ 
ing values. However, we strongly prefer models with g 
= 0: those with g > 0.05 predict disk ansae that are 
too faint compared to regions at small angular separa¬ 
tions. Models with e < 0.1 or e = 0.3 fail to prop¬ 
erly trace the observed disk. A single model with e 
= 0.1 fits our threshold, while all other acceptedly- 
fitting models have e=0.2, suggesting a disk eccentric¬ 
ity (e = 0.1-0.2) compar a,ble to that of the most eccen¬ 
tric k nown debris disks (|Kalas et al.l 120051 : IKrist et all 
1201211 . The deprojected, deconvolved, and azimuthally- 
averaged best-fit disk model has a normalized FWHM of 
0.37, among the largest for scattered light-imaged debris 
rings (|Rodigas et al.l[2Mll . 

4.3. Dynamical Sculpting of the Disk: Limits on Planets 

Debris disk dynamical modeling sets stronger limits 
on unseen planets than from the GPI contrasts limits 
alone (contr ast > 10~^ at r < 0 74 (44 AU) or Mpianet 
> 7 Mj; c.f. iBaraffe et al.l 1200311 . The eccentric ring is 
plausibly due to dynamical perturbations from an un¬ 
seen planet, which can open a gap in the disk and stir 
it via secular perturbations. In both cases, the planet’s 
eccentricity can be estimated assuming that the eccen¬ 
tricity of th e ring equals the local value of the forced 
eccentricity (IQuillenl 1200611 . First, we place limits on 
the properties of the planet sculpting the disk from the 
iNesvold and Kuchnerl (|2015H gap opening model, where 
the planetesimals are initially dynamically hot (0.0 < e 
< 0.2) and the planet interacts with the disk via mean- 
motion resonances, triggering collisions and producing an 
inner hole. Gap opening depends on the system age in 
units of the grain collisional time, which in turn depends 
on the orbital period and op tical depth: fnrhit,/(4 7r r). 

Second, we explore the iMustill and Wvatti (|2009(1 
planet stirring model: the planet stirs an initially cold 
disk from the inside-out via secular interactions, and the 
radius of the central hole extends to where planetesimal 
orbits begin to cross. We compute stirring from first- 
order secular theory (valid for low planet eccentricities), 
corresponding here to larger planet semimajor axes. Our 
two calculations assume 48 AU as the ring center, the 
fractional luminosity as a proxy for optical depth, and 
an age of 15 Myr. 

Interior to ^ 30 AU only a superjovian-mass planet can 
sculpt the ring by gap opening (Figure [2|). Planets with 
masses and semimajor axes comparable to the outer solar 
system planets could stir the disk to appear as a bright 
debris ring with an eccentricity of 0.1-0.2 (e.g. a Saturn 


with e = 0.2). In both cases, Super-Earths just interior 
to the ring edge could sculpt the ring. 

4.4. Spectral Analysis 

We perform two spectral extractions. First, we ex¬ 
tract the raw spectrum at both disk ansae, identifying 
the ‘centroid positions from the collapsed image. We 
then extract the spectrum of the best-fit disk model at 
these centroid positions. In both cases, we measure the 
surface brightness: the mean brightness within the same 
8-pixel diameter aperture used for flux calibration with 
errors calculated using method 2 described in Section 3. 
For the model disk, we average the results from the two 
ansa and correct for signal loss due to a finite aperture 
by comparing the surface brightness of the convolved and 
unconvolved best-fit disk models. 

Figure [3^ shows our extracted best-fit scaled model 
disk spectrum (blue), averaged raw spectrum (green), 
and raw spectra extracted for each ansa (dotted lines). 
The spectra have a relatively flat to slightly blue slope 
across H band. Spectra extracted from the two ansae 
exhibit strong agreement. 

The debris disk’s scattered light spectral shape is sen¬ 
sitive to its dust’s co mposition. Inste ad of using ad¬ 
vanced methods (e.g. iLisse et HI I2007D . we model the 
reflectance of dust using simple Mie theory-based sin¬ 
gle constituents found amongst dust and large bodies 
in the Kuiper belt- carbonaceous d ust, amorphous si l- 
icates, and water ice - drawn from ILisse et al.l (II998I1 . 
nominally using the best-fit Halley-like particle size dis¬ 
tribu tion ([Krishna Swamv and ShahllI988t iMazets et al.l 
Il986f) but exploring simple power laws with exponents 
between -3 and -4. To calculate reflectance, we divide 
the best-fit model disk’s spectrum by that of a Pickles 
F2V star. 

The raw (black curve) and binned (cyan curve) disk 
reflectance spectra appear neutral/slightly blue at 1.5- 
1.75 /im (Figure [Sja). Reflectance spectra of water ice 
(blue), amorphous silicates (green) and amorphous car¬ 
bon (red) show differences across the H passband. Wa¬ 
ter ice (and silicates) with a Halley-like size distribution 
provide the best match, reproducing the spectrum’s neu¬ 
tral/blue slope. Amorphous carbon has a x^ value 2.5 
times higher than water ice’s and is marginally disfa¬ 
vored (at the 68% confidence limit) for the binned spec¬ 
trum, having a goodness-of-fit statistic of 0.78. However, 
adopting a simple power law exponent for the particle 
size distribution between —3 and —4, a subset of mod¬ 
els from each species likewise fit. Thus, our spectrum’s 
large error bars and narrow wavelength range currently 
preclude us from making definitive statements about the 
dust composition. 

Comparing the H-hand spectrum with thermal IR data 
and properties of other debris rings may better reveal 
evidence for ice. HD 115600’s thermal disk emission 
(Figure Ej:) peaks at 115 K, but blackbody-like emis¬ 
sion should be located at 13 AU, not the observed 48 
AU. Real dust around stars later than AO should be 
located at increasingly larger distances than predicted 
by blackbody emission , over 2.5 times fur ther for an F2 
star (see Figure 10 in iBooth et al. 11201^ or beyond 33 
AU in HD 115600’s case. Steeper size distributions and 
sub-blowout dust sizes (expected for a collisionally-active 
disk) may be sufficient to yield dust around HD 115600 
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at the right distances. However, organics (carbon) dom¬ 
inated dust yields dust locations an additional factor of 
4-8 larger than predicted by blackbody emission; ice and 
ice/silicat e mixtures yield ne gligible enhancement (Fig¬ 
ure 12 in iBooth et al. I[2ni3[ l. A water-ice (dominating 
scattered light) and organics (dominating thermal emis¬ 
sion) mixture, much like that found on Kuiper belt object 
surfaces, may fit as well. 

HR 4796A’s f ractional luminosity is over 2 times higher 
(- 4.8 X 1 n-3: IT.ow et~ZI[2nn^ IChen et abll^nTl . The 
i7-band surface brightness contrast of HD 115600’s disk 
ansae is ^ 1.6 mag nitudes brighter than HR 4796A’s 
(|Rodigas et al.l[201^ . Moving the HR 4796A ring from 
70 AU to HD 115600’s disk location (48 AU) still re¬ 
sults in scattered light emission about half as bright as 
HD 115600’s disk. Thus, HD 115600’s disk is reflecting 
light more efficiently than HR 4796A’s disk while having 
less thermal emission, a result explicable if HD 115600s 
disk is dominated by higher albedo species like water-ice. 
Multi-wavelength photometry/spectroscopy is needed to 
more decisively assess the composition of HD 115600’s 
disk. 

5. DISCUSSION 

HD 115600’s debris ring is the first object dis¬ 
covered using the new generation of extreme adap- 
tive optics instruments (SCExAO, GPI, and SP HERE; 
Martinache and GuvonI 120091 : [Macintosh et al.l 12014 
Beuzit et al.l l2008l) and the first debris ring with 
spatially-resolved integral field spectroscopy. The ring 


is confined to a Kuiper belt-like stellocentric distance 
and its reflectance is consistent with that of major con¬ 
stituents of Kuiper belt bodies and their ejecta, including 
water ice. As the disk orbits a 15 Myr-old star only ^ 
40-50% more massive than the Sun located in the nearest 
OB association, it provides a promising reference point 
for understanding the early evolution and composition of 
the Kuiper belt. 

HD 115600’s disk shows evidence for dynamical sculpt¬ 
ing by a solar system-like giaiit plan et. With the excep¬ 
tion of P Pic b (|Lagrange et al.l[201fl[ l. fully-formed young 
planets directly imaged thus far have super-jovian masses 
(5-10 Mj) and orbit at wide separations (15-150 AU) 


(iMarois et al.ll2008|: 

Rameau et al.ll2013l:lKuzuhara et al.l 

120131: ICurrie et al.f 

2014ar). Southern hemisphere svs- 


terns GPI and SPHERE could image even lower-mass 
planets located just interior to HD 115600’s debris ring 
and superjovian planets near their inner working an¬ 
gles (r ^ 0'.'l-0'.'2 or ^ 10-20 AU). Combining new 
planet detections/upper limits with more precise esti¬ 
mates of the disk’s properties will make HD 115600 an 
excellent laboratory for studying planet- disk interactions 
(|Chiang et al.ll200a iRodigas et al.l[2?)T^ . 
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TABLE 1 

Debris Disk Forward Modeling 


Parameter 

Model Range 

Best-Fit Model 

Well-Fitting Models 

ro (AU) 

47 . 

.. 48 

48 

47 . 

.. 48 

e 

0 .. 

. 0.3 

0.2 

0.1 

... 0.2 


7.5 

... 10 

7.5 

7.5 

... 10 

(^out 

-5 .. 

. -7.5 

-7.5 

-5 

... -7.5 

Ay (AU) 

0 . 

.. 1 

1 

0 ... 

1 

9 

0 ... 

0.15 

0 

0 


ksio (AU) 

1., 

.. 5 

3 

1 ... 

5 
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Fig. 1.— (top-left) Reduced, PSF-subtracted, and wavelength-collapsed image of HD 115600 showing a bright, nearly edge-on debris ring 
located just beyond a Pluto-like distance to the star. The disk center (diamond), determined from ellipse fitting, appears offset from the 
star’s position (cross), (top-right) A representative wavelength slice (1.671 fim) of the final datacube. (bottom-left) The best-fit model disk 
produced by GRaTeR: ao = 48 AU, PA = 24°, i = 79.5°, Oo = 90°, ksio = 3 AU, g = 0, e = 0.2, ain = 7.5, and aout = —5. (bottom-right) 
The same model, convolved with the appropriate PSF in each channel, PSF subtracted, and wavelength-collapsed. The model reproduces 
the suppressed signal in the northeast and southwest portions of the disk, although there are some small discrepancies (see main text). 
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Fig. 2.— Limits on the mass of a planet stirring ring planetesimals or opening a gap in the disk and clearing material out into a ring. 
Planets much like the outer solar system planets (overplotted) can stir a disk with an eccentricity consistent with the ring properties. 
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Fig. 3. — (Top-Left) Surface brightness of the best-fit model disk spectrum (averaged between ansa). Due to low throughput, we trim 
the first and last two channels. The raw disk spectrum shows an almost identical shape but consistently 15% lower signal: spectra 
extracted from individual ansae (dashed curves) agree within their error bars. (Top-Right) Reflectance spectra of the best-fit model disk 
and a binned (to the spectral resolution of GPI) version of the best-fit mod el compared to Mi e theory predictions for water ice, amorphous 
silicates, and amorphous carbon, (bottom) HD 115600’s optical/IR SED IjChen et al.ll2Q15h revealing 115 K ± 6 K (2-c7) cold dust and 
favoring an icy/ice-silicate composition (see main text). 



















































